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Abstract. The final results of the MACRO experiment on atmospheric neutrino oscillations are presented. 
The data concern difTerent event topologies with average neutrino energies of ~ 3 and ~ 50 GeV. Multiple 
Coulomb Scattering of the high energy muons was used to estimate the neutrino energy event by event. The 
angular distributions, the L/E^ distribution, the particle ratios and the absolute fluxes all favour Vf^ -^ Vr 
oscillations with maximal mixing and Am'' ~ 0.0023 eV^. Emphasis is given to measured ratios which are 
not affected by Monte Carlo (MC) absolute normalization; a discussion is made on MC uncertainties. 
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1 Introduction 

MACRO was a large area multipurpose underground de- 
tector ^ designed to search for rare events and rare phe- 
nomena in the penetrating cosmic radiation. It was lo- 
cated in Hall B of the underground Gran Sasso Lab at an 
average rock overburden of 3700 m.w.e.; it started data 
^ taking with part of the apparatus in 1989; it was com- 
Si^ ■ pleted in 1995 and was running in its final configuration 
5— I ' until the end of 2000. The detector had global dimensions 
^ of 76.6 X 12 X 9.3 m'^ and provided a total acceptance to 
an isotropic flux of particles of ^ 10, 000 vc? sr; vertically 
it was divided into a lower part, which contained 10 hor- 
izontal layers of streamer tubes, 7 of rock absorbers and 

2 layers of liquid scintillators, and an upper part which 
contained the electronics and was covered by 1 layer of 
scintillators and 4 layers of streamer tubes. The sides were 
covered with 1 vertical layer of scintillators and 6 of lim- 
ited streamer tubes. 

MACRO detected upgoing I'^'s via charged current inter- 
actions, v^ — > [i\ upgoing muons were identified with the 
streamer tube system (for tracking) and the scintillator 
system (for time-of- flight measurement). The events mea- 
sured and expected for the three measured topologies, de- 
viate from Monte Carlo expectations without oscillations, 
FigE these deviations and the L/Ey distribution point to 
the same v^ — > Vr oscillation scenario |2]-[H|, Fig. |21 



2 Atmospheric neutrinos. Monte Carlo 

The measured data of Fig. ^were compared with differ- 
ent MC simulations. In the past we used the neutrino flux 



computed by the Bartol96 group ^ and the GRV94 par- 
ton distribution. For the low energy channels the cross 
sections by P. Lipari et al. were used; the propagation 
of muons to the detector used the energy loss calculation 
by Lohmann et al. The total systematic uncertainty in the 
predicted flux of upthroughgoing muons, was estimated at 
±17 %; this is mainly a scale error that does not change 
the shape of the angular distribution. The response of the 
detector to the passage of particles was simulated using 
GEANT3. A similar MC (Honda96) was used by the Su- 
perK Collaboration [E] [TI] . 

Recently new improved MC predictions for neutrino fluxes 
were made available by the Honda ^J and FLUKA jT2] 
groups. They include three dimensional calculations of 
hadron production and decays and of neutrino interac- 
tions, improved hadronic model and new fits of the pri- 
mary cosmic ray flux. The two MC yield predictions for 
the non oscillated and oscillated v^ fluxes equal to within 
few % ^, . The shapes of the angular distributions for os- 
cillated and non oscillated Bartol96, new FLUKA and new 
Honda fluxes are the same to within few %. The absolute 
values of our upthroughgoing muon data are about 25% 
above those predicted by the new FLUKA and Honda MC, 
Fig. El A similar situation is found in the new SuperK data 
[inj . The high energy i/^ data thus suggest that the new 
Honda and FLUKA predictions are low, probably because 
of the used CR fit (the inclusion of the new ATIC CoUab. 
measurements of primary CRs may improve the situation 
|13|). The evidence for neutrino oscillations is due mainly 
to the shape of the angular distribution and this is the 
same in all MCs. Also the ratios of the medium to high 
energy measurements and of the two different samples of 
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low energy data are MC independent. Our low energy data 
suggest that the FLUKA normalization should be raised 
by about 12% at E^ r^ 3 GeV. 
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Fig. 1. Zenith distributions for the MACRO data (black 
points) for (a) upthroughgoing, (b) semicontained and (c) up- 
stopping muons + down semicontained. The dashed line in (a) 
is the no-oscillation Bartol96 MC prediction (with a scale error 
band), (in b, c are the FLUKA MC; better fits are obtained 
with Bartol96); the solid lines refers to i/^ — > Vt oscillations 
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Fig. 2. Ratio (Data/MC Bartol96) versus the estimated L/Ei, 
for the upthrougoing muon sample (black circles) and the semi- 
contained up-^ (open circle). The horizontal dashed line at 
Data/MC=l is the expectation for no oscillations. 
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Fig. 3. Comparison of our measurements with the Bartol96 
and the new Honda 2001 oscillated and non oscillated fluxes. 



the whole detector. The corresponding f^/s have a median 
energy Ei, ~ 50 GeV. Many possible systematic effects 
and backgrounds that could affect the measurements were 
studied |31|S]. The data, Fig. ^, deviate in shape and in 
absolute value from the Bartol96 MC non oscillated pre- 
dictions. 

j/^ -^ v-r versus Vf^ —^Vg- Matter effects would produce a 
different total number and a different zenith angle distri- 
bution of upthroughgoing muons. The ratio Ri — Ver- 
tical/Horizontal = iV(-l < cos9 < -0.7)/N{-0A < 



with maximal mixing and Am = 2.3 ■ 10" eV (see text). ^^^q ^ q) ^^^ ^^^^ ^^ ^^^^ ^^^ ^^ ^ ^^ oscillation hypoth 



esis versus v,. 



0. The 



3 MACRO results on atmospheric neutrinos 

The upthroughgoing muons come from f^ interactions in 
the rock below the detector; muons with _E^ > 1 GeV cross 



oscillations 

(with any mixing) are excluded at about 99.8% c.l. with 
respect to v^ — > Vt oscillations with maximal mixing |H]. 
Oscillation probability as a function of the ratio L/E^. 
Ejj was estimated by measuring the muon energy, _E^, by 
means of the muon Multiple Coulomb Scattering (MCS) in 
the rock absorbers in the lower MACRO. The space reso- 
lution achieved is ~ 3 mm. For each muon, seven variables 
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were given in input to a Neural Network (NN) trained to 
estimate muon energies with MC events of known input 
energy crossing the detector at different zenith angles. The 
distribution of the ratio R = [Data/MCnoosc) obtained by 
this analysis is plotted in Fig. |21as a function of {L/E„) 
13. The data extend from {L/E^) - 30 to 5000 km/GeV. 




Fig. 4. Interpolated qualitative 90% C.L. contour plots of the 
allowed regions in the Am^ — sin^ 26 plane for the MACRO 
data using only the ratios Ri,R2, R3 (outer continuous line) 
and using also the absolute values assuming the validity of the 
Bartol96 fluxes (dotted line). 



The Internal Upgoing (lU) muons come from ~ 3 GeV 
i^^'s interacting in the lower apparatus. Compared to the 
no-oscillation prediction there is a reduction in the flux of 
these events, without distortion in the shape of the zenith 
distribution, Fig. ^p. The MC predictions for no oscilla- 
tions in Figs.n^ and^ are the dashed lines with a 21 % 
systematic band. At these energies the ratios DATA/MC 
are 0.67 for FLUKA and 0.54 for Bartol96. 

The Upstopping (UGS) muons are due to ~ 3 GeV I'^'s 
interacting below the detector, yielding upgoing muons 
stopping in the detector. The Semicontained Downgoing 
(ID) muons are due to I'^-induced downgoing /i's with 
vertex in the lower MACRO. The two types of events are 
identified by topological criteria. The upgoing i/^'s should 
have oscillated completely, while the downgoing v^ do not. 



4 Determination of the oscillation 
parameters. 

In the past, in order to determine the oscillation param- 
eters, we made fits to the shape of the upthroughgoing 
muon zenith distribution and to the absolute flux com- 
pared to the Bartol96 MC prediction. The other data were 
only used to verify the consistency and to make checks. 
The result was Am? = 0.0025 eV^ and maximal mixing 
in] [3. Later, also the L/E^ distribution was considered 

m 

In order to reduce the effects of possible systematic 
uncertainties in the MC we recently used the following 



three independent ratios [H| and we checked that FLUKA, 
Honda and Bartol96 Monte Carlo simulations yield the 
same predictions to within ^ 5%. 

(i) High Energy Data: zenith distribution ratio: i?i = 

Nvert/Nhor 

(ii) High Energy Data, neutrino energy measurement ra- 
tio: i?2 = Nlow/Nhtgh 

(iii) Low Energy Data: 

Ratio i?3 = {Data/MC)iu/{Data/MC)iD+UGS- 

With these ratios, the no oscillation hypothesis has a prob- 
ability P ^ 3 ■ 10~^ and is thus ruled out by ^ 5(j. By 
fitting the three ratios to the v^ —>■ v-r oscillation formu- 
lae we obtain sin^ Id = I, Arr? = 2.3 • lO'^ eV^ and the 
allowed region indicated by the solid line in Fig. 0] 

There is a good consistency between the old and new 
methods. 

If we use the Bartol96 flux we may add to the ratios 
(i) - (iii) the information on the absolute flux values of the 

(iv) high energy data (systematic scale error of ;^17%) i?4 = 

Nmeas /NMCBartol ■ 

(v) low energy semicontained muons, with a systematic 

scale error of 21%, Rr, ~ Nmeas/NMCBartol- 

These informations reduce the area of the allowed region 
in the Am^ — sin^ 29 plane, as indicated by the dashed 
line in Fig. ^ The limit lines represent smoothed inter- 
polations and are qualitative. The final MACRO Am'^ is 
2.3-10-3 eV2. 

We would like to acknowledge the cooperation of the mem- 
bers of the MACRO collaboration. 
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